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Abstract 

t/3 , We present an analytical framework to assess the link layer throughput of multichannel Opportunistic 

O 

Spectrum Access (OSA) ad hoc networks. Specifically, we focus on analyzing various combinations of 

CO ' collaborative spectrum sensing and Medium Access Control (MAC) protocol abstractions. We decompose 

> 

' collaborative spectrum sensing into layers, parametrize each layer, classify existing solutions, and pro- 

O ' 

, pose a new protocol called Truncated Time Division Multiple Access (TTDMA) that supports efficient 

distribution of sensing results in "k out of N" fusion rule. In case of multichannel MAC protocols we 
evaluate two main approaches of control channel design with (i) dedicated and (ii) hopping channel. We 
propose to augment these protocols with options of handling secondary user (SU) connections preempted 
by primary user (PU) by (i) connection buffering until PU departure and (ii) connection switching to 



o 



. a vacant PU channel. By comparing and optimizing different design combinations we show that (i) it 

' is generally better to buffer preempted SU connections than to switch them to PU vacant channels and 

(ii) TTDMA is a promising design option for collaborative spectrum sensing process when k does not 
change over time. 
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I. Introduction 

It is believed that Opportunistic Spectrum Access (OSA) networks will be one of the primary forces in 
combating spectrum scarcity lH in the upcoming years 131, H. Therefore, OSA networks Q, ||6l have 
become the topic of rigorous investigation by the communications theory community. Specifically, the 
assessment of spectrum sensing overhead on OSA medium access control (MAC) performance recently 
gained a significant attention. 

A. Research Objective 

In the OSA network performance analysis, a description of the relation between the primary (spectrum) 
user (PU) network and the secondary (spectrum) user (SU) network can be split into two general models: 
macroscopic and microscopic. In the macroscopic OSA model HI, S it is assumed that the time 
limit to detect a PU and vacate its channel is very long compared to the SU time slot, frame or packet 
length duration. Such a time limit is assumed to be given by a radio spectrum regulatory organization. 
For example, the timing requirements for signal detection of TV transmissions and low power licensed 
devices operating in TV bands by IEEE 802.22 networks iflOl (including transmission termination and 
channel vacancy time, i.e. a time it takes the SU to stop transmitting from the moment of detecting PU) 
must be equal to or smaller than 4.1 s ifTTl Tab. 15.5], while the frame and superframe duration of IEEE 
802.22 are equal to 10ms and 160ms, respectively iTm . Also, in the macroscopic model it is assumed 
that the PU channel holding time, i.e. the time in which the PU is seen by the SU as actively transmitting, 
is much longer than the delay incurred by the detection process performed at the SU. As a result it can 
be assumed in the analysis that, given high PU detection accuracy (which is a necessity), OSA network 
performance is determined by the traffic pattern of the SUs. That is, it depends on the total amount of 
data to be transmitted by the SU network, the duration of individual SU data packets and the number 
of SU nodes. In other words the PU bandwidth resource utilization by the SU is independent of PU 
detection efficiency. 

In the microscopic OSA model, more popular than its macroscopic counterpart due to analytic chal- 
lenges, the detection time is short in relation to the shortest transmission unit of the OSA system. 
Detection is also performed much more frequently than in the macroscopic model, i.e. for every SU 
packet |[T2l . |[T3l or in every time slot |[T4l . (15\ . |[T6ll . IfTTl . ifTSl . Also, the microscopic model assumes 
much higher PU activity than the macroscopic model, which justifies frequent detection cycles. Since the 
detection overhead is much larger than in the macroscopic model, the analysis of utilization of resources 
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(temporarily unoccupied by PU) by OSA network cannot be decoupled from the analysis of the PU signal 
detection phase. 

Therefore, while the distinction between macroscopic and microscopic models are somehow fluid, it 
is important to partition the two cases and compare them in a systematic manner. More importantly, the 
comparison should be based on a detailed OSA multichannel and multiuser ad hoc network model |fT9l 
Sec. 7.4], which would not ignore the overhead from both the physical layer (PHY) and MAC layers of 
different cooperative and distributed spectrum sensing strategies |[T9l Tab. 7.1] and, in case of microscopic 
model, account for different channel access procedures and connection management strategies for the SUs 
upon PU detection, like buffering or switching to a vacant channel. Finally, the comparison should be 
realized using tractable analytical tools. 

B. Related Work 

The literature on this topic can categorized into three groups: (i) performance analysis of general 
OSA networks, excluding a detailed model for spectrum sensing (mostly for the macroscopic model), 
(ii) performance of spectrum sensing isolated from MAC aspects of network collaboration, and (iii) joint 
performance of spectrum sensing and networking for OSA. 

One of the first works that gained insight into the general performance of OSA networks, considering 
impact of PU activity on blocking and throughput of the SU network was |[20l . where the capacity 
of a multichannel OSA system was assessed by comparing centrally coordinated versus random SU 
channel assignment. A spectrum sensing process was not considered. A similar problem was investigated 
in |[2n where the spectrum sharing gains for PU and SU networks were obtained for a distributed and 
multichannel ad hoc OSA network. Unfortunately, a zero delay spectrum sensing process was assumed 
with genie-aided channel selection, i.e. in every time slot the receiver knew of the exact channel the 
transmitter will use to send data ll2ll, Sec. III-Cl]. 

In later works, assumptions on the OSA network model became more realistic. Specifically, Markovian 
analysis of SU traffic buffering on the event of PU arrival was presented for a SU exponential service 
time |22| and for a SU phase-type service time |[23l . Unfortunately the impact of spectrum sensing 
detection time overhead on the OSA network performance was not investigated and the connection 
arrangement process for new SU arrivals, i.e. method to select and access a channel for a new sender- 
receiver pair, was assumed to be performed by a centraUzed entity. A different option of the above model 
has been analyzed in |24|, with only PU channels dedicated to the OSA network and with a mixture 
of PU and SU exclusive channels. SU connection buffering was not allowed, however, SU connections 
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were able to switch to an empty SU exclusive channel on the event of channel preemption by the PU. 

A similar analysis, but with a different channelization structure, where the PU occupied more than 
one SU channel (contrary to |[22l . |[23l ) was performed in l25l . The authors addressed the cases of (i) 
connection blocking, and (ii) channel reservation and switching of SU connections to empty channels on 
PU arrival. This analysis was later extended to the case of finite SU population and packet queuing |26|, 
and buffering and switching of SU connections preempted by PU arrivals |27|. Again, in all papers listed 
above the spectrum sensing process was assumed to have no overhead and perfect reliability. Moreover 
the connection arrangement process for SUs was not considered. 

A system where the PU had to wait until a SU vacates a channel was analyzed in |[28l . Both perfect and 
imperfect PU detection processes were considered, however detection overhead as well as a connection 
arrangement process for the secondary system was not considered. In 1291 another OSA system consisting 
of only two SUs was considered. Once again, SU connection arrangement process was not included in 
the analysis and two distinct strategies were assumed: (i) with complete SU connection blocking on the 
arrival of PU (just like in Il24l . |[25l ) and (ii) with SU connection buffering on the arrival of PU (just 
like in |[27l ). Spectrum sensing and connection arrangement process was not investigated. A practical 
OSA model analyzing the impact of PU activities on the quality of Voice over Internet Protocol traffic 
was analyzed in |[30l with a system not allowing SU connection buffering and perfect information on 
the PU channel presence assumed. The only work considering a microscopic model was |[T5l . where 
a relation between sensing time, PU detectability and different connection arrangement processes were 
considered. Detailed simulations and Markov analysis were performed, however, as noted in the paper, 
the proposed model did not yield accurate results over the ranges of all parameters considered, e.g. level 
of PU activity. Also, only one sensing strategy with SU connection buffering was analyzed for the case 
of different MAC protocols. 

Considering a second group of papers (related to the performance of spectrum sensing algorithms in 
isolation from higher protocol layers), in llT4l the analysis of the average time consumed by two stage 
spectrum sensing (proposed independently in ||8l, lITSl . |[3T1 ) was decoupled from SU traffic characteristics. 
Moreover the delay caused by exchanging hard decision measurements in the cooperative sensing process 
based on the "and" rule was not included. The impact of sensing overhead on, e.g., throughput, and energy 
consumption was explored in |[32l . However the analysis did not account for any OSA network and SU 
traffic. Also, the relationship between detection time and detection quality was not investigated. In fTSl 
a microscopic model was analyzed with a sensing period every slot, synchronization between PU and 
SU, and PU stationary over the whole slot duration. Markov analysis was performed only to evaluate 
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the delay incurred while searching for unoccupied spectrum. The sensing process was not coupled to 
any of the known MAC protocols and SU connections. Also, only non-collaborative spectrum sensing 
was considered. Finally, in fj\ most of the procedures related to spectrum sensing were categorized and 
divided into Open Systems Interconnection-like layers. Performance of the most common combinations 
of sensing algorithms were assessed, but only for the macroscopic model. Our analysis is the microscopic 
treatment of Q. 

Considering the final group of papers, when coupling spectrum sensing procedures with link layer 
protocols, there is a fundamental tradeoff between sensing time, sensing quality and OSA network 
throughput. This has been independently found for general OSA network models with a single sensing 
band |[T2l . multiple sensing bands ll33l with and without cooperative detection and centralized resource 
allocation, and in a context of MAC protocol abstraction |[T5l for a non-cooperative sensing case. See also 
recent discussion in |[T9l Sec. 2.3.1, 7.3, and 10.2.4]. This tradeoff will be especially clear, while evaluating 
microscopic models, since the detection time creates a significant overhead for the data exchange phase. 
Recently the model of |12| was extended to the case of "k out of N" rule in cooperative sensing [131, 
optimizing parameters of the model to maximize the throughput given detection rate requirements. 
Unfortunately, the delay caused by exchanging sensing information was not included. 

C. Our Contribution 

In this paper, we present a unified analytical framework to design the spectrum sensing and the OSA 
data MAC jointly, for the macroscopic and microscopic cases. This design framework provides the (i) 
means of comparing different spectrum sensing techniques plus MAC architectures for OSA networks and 
(ii) spectrum sensing parameters such as observation time and detection rate for given design options. As 
a metric for optimization and comparison, we consider the average link layer OSA network throughput. 
Our model will account for the combined effects of the cooperative spectrum sensing and the underlying 
MAC protocol. For spectrum sensing, we will consider several architectures parametrized by sensing 
radio bandwidth, the parameters of the sensing PHY, and the parameters of the sensing MAC needed to 
exchange sensing data between individual OSA nodes. Along with classifying most of the well known 
sensing MAC protocols, we introduce a novel protocol called Truncated Time Division Multiple Access 
(TTDMA) that supports efficient exchange of individual sensing decisions in "k out of A^" fusion rule. 
For the data MAC we will consider two protocol abstractions, (i) Dedicated Control Channel (DCC) and 
(ii) Hopping Control Channel (HCC), as analyzed in |[T5l . |[34|| with novel extensions. That is, given 
the designs of 1251, 1261, 123, l30l, we will analyze MAC protocols that (i) allow (or forbid) to buffer 
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existing SU connections on the event of PU arrival, and (ii) allow (or forbid) to switch the SU connections 
preempted by the PU to the empty channels. Please note that in the case of the analytical model proposed 
in |[T5l for the SU connection buffering OSA MAC schemes we present an exact solution. Finally, using 
our framework, we compute the maximum link layer throughput for most relevant combinations of 
spectrum sensing and MAC, optimizing parameters of the model jointly, both for the microscopic and 
macroscopic models. 

The rest of the paper is organized as follows. System model and a formal problem description is 
presented in Section JIl Description of spectrum sensing techniques and their analysis is presented in 
Section |llll Analysis of MAC strategies are presented in Section |lVl Numerical results for spectrum 
sensing process, MAC and joint design framework are presented in Section |Vl Finally the conclusions 
are presented in Section |Vll 

II. System Model and Formal Problem Description 

The aim of this work is to analyze link layer throughput accounting for different combinations of 
MAC, spectrum sensing protocols and regulatory constraints. The model can later be used to optimize 
the network parameters jointly to maximize the throughput, subject to regulatory constraints. Before 
formalizing the problem, we need to introduce the system model, distinguishing between the microscopic 
and macroscopic approaches. 

A. System Model 

1) Microscopic Model: For two multichannel MAC abstractions considered, i.e. DCC and HCC, we 
distinguish between the following cases: (i) when SU data transfer interrupted by the PU is being buffered 
(or not) for further transmission and (ii) when existing SU connection can switch (or not) to a free channel 
on the event of PU arrival (both for buffering and non-buffering SU connection cases). Finally, we will 
distinguish two cases for DCC where (i) there is a separate control channel not used by the PU and (ii) 
when control channel is also used by the PU for communication. All these protocols will be explained 
in detail in Section ITVl 

We assume slotted transmission within the SU and PU networks, where PU and SU time slots are equal 
and synchronized with each other The assumptions on slotted and synchronous transmission between PU 
and SU are commonly made in the literature, either while analyzing theoretical aspects of OSA (see lfT2l 
Fig. 2], lUSl Sec. Ill], im Fig. 1], (SU Fig. 1], (Ml Fig. 5 and Sec. 5.2], gTl Fig. 1]) or exploring practical 
OSA scenarios (see |[T6l Fig. 2] in the context of secondary utilization of GSM spectrum or |[38l in the 
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context of secondary IEEE 802.16 resources usage). Our model can be generalized to the case where PU 
slots are offset in time from SU slots, however, it would require additional analysis of optimal channel 
access policies, see for example |[36l . |[39]| . BOl . which is beyond the scope of this paper. We also note 
that the synchrony assumption allows one to obtain upper bounds on the throughput when transmitting 
on a slot-asynchronous interface BTI . 

The total slot duration is tt /US. It is divided in three parts: (i) the detection part of length tq ^us, denoted 
as quiet time, (ii) the data part of length t„ /us, and if communication protocol requires channel switching 
(iii) switching part of length tp /is. The data part of the SU time slot is long enough to execute one 
request to send and clear to send exchange |[T5l . ll34l . For the PU the entire slot of tt fis is used for data 



transfer, see Fig. 1(a) 



Our model assumes that there are M channels having fixed capacity CMbps that are randomly and 
independently occupied by the PU in each slot with probability qp. There are N nodes in the SU network, 
each one communicating directly with another SU on one of the available PU channels in one hop fashion. 
Also, we assume no merging of the channels, i.e only one channel can be used by a communicating pair of 
SUs at a time. SUs send packets with geometrically distributed length with an average of 1/g = d/{Ctu) 
slots for DCC, and l/q = d/{C{tu + tp}) slots for HCC US] Sec III-C4b], iMl Sec. 3.2.3], where d is 
the average packet size given in bits. Difference between average packet length for DCC and HCC is a 
result of switching time overhead for HCC, because during channel switching SUs do not transfer any 
data, even though they occupy the channel. We therefore virtually prolong data packet by tp for HCC to 
keep the comparison fair. 

Every time a node tries to communicate with another node it accesses the control channel and transmits 
a control packet with probability p to a randomly selected and non-occupied receiver. A connection is 
successful when only one node transmits a control packet in a particular time slot. The reason for selecting 
a variant of S -ALOHA as a contention resolution strategy was manyfold. First, in reality each real-life 
OSA multichannel MAC protocol belonging to each of the considered classes, i.e. HCC or DCC, will 
use its own contention resolution strategy. Implementing each and every approach in our analysis: (i) 
would complicate significantly the analysis, and most importantly (ii) would jeopardize the fairness of 
the comparison. Therefore a single protocol was needed for the analytical model. Since S-ALOHA is a 
widespread and well understood protocol in wireless networks and is a foundation of many other collision 
resolution strategies, including CSMA/CA, it has been selected for the system model herein. 

In each quiet phase every SU node performs PU signal detection based on signal energy observation. 
Since we assume that OSA nodes are fully connected in a one hop network, thus each node observes on 
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average the same signal realization in each time slot llT3l . lITSl . ll42l . PU channels detected by the SU are 
assumed as Additive White Gaussian Noise with a channel experiencing Rayleigh fading. Therefore to 
increase the PU detectability by the OSA network we consider collaborative detection with hard decision 
combining in the detection process based on "k out of N" rule, as in ll43l . fl4l . Hence we divide the 
quiet phase into (i) the sensing phase of length tg /xs and (ii) the reporting phase of length tr /iS. The 
sensing phase is of the same length for all nodes. For simplicity we do not consider in this study sensing 
methods that adapt the sensing time to propagation conditions as in ll45l . In the sensing phase, nodes 
perform their local measurements. Then, during the reporting phase, nodes exchange their sensing results 
and make a decision individually by combining individual sensing results. We will analyze different 
PHY and MAC approaches to collaborative spectrum sensing, especially (i) methods to assign sensing 
frequencies to users, (ii) rules in combining the sensing results, and (iii) multiple access schemes for 
measurement reporting. In this paper we do not consider sensing strategies applicable to single channel 
OSA networks H6l . two stage spectrum sensing El, and sensing MAC protocols based on random 
access pTl . due to their excessive delay. We will explain our spectrum sensing approaches in more detail 
in Section |llll Further we assume a error channel, for the sensing layer as well as for data layer where 
probability of error during transmission is denoted as p^. 

Finally, we consider two regulatory constraints under which the OSA network is allowed to utilize 
the PU spectrum provided the channel is idle: (i) maximum detection delay td,max, i e. a time limit 
within which a SU must detect a PU, and (ii) minimum detection probability pd,mm, i-c, a probability 
with which a OSA system has to detect a PU signal with minimum signal to noise ratio 7. Note that 
in the event of mis-detection and subsequent SU transmission in a channel occupied by PU, a packet 
fragment is considered successfully transmitted, since in our model transmission power of SU is much 
higher than interference from PU, and regulatory requirements considered here do not constrain SU 
transmission poweiQ (refer for example to IEEE 802.22 draft where Urgent Coexistent Situation packets 
are transmitted on the same channel as active PU ifTOl . ifTTI ). Moreover, maximum transmission power is 
a metric specific to overlay OSA systems |fT9l Sec. 2.2.5 and 8.2.1] where typically no spectrum sensing 
is considered. Also we do not consider a metric based on a maximum allowable level of collisions 
between PU and SU. Note that the parameters of the introduced model are summarized in Table U and 

'The opposite case is to assume that a packet fragment is considered as lost and retransmitted. This approach however requires 
an acknowledgement mechanism for a lost packet fragment, see for example 1171 Sec. II-B3], |41 Sec. II], that contradicts the 
model assumption on the geometric distribution of SU packets. 
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Fig. 1. Difference between macroscopic and microscopic model in relation to sensing process and slot structure: (a) microscopic 
case and (b) macroscopic case. Symbols are explained in Section Hl-AI and Table |I] 



the abbreviations are summarized in Table |lll 

2) Macroscopic Model: We assume the same system model as for the microscopic case, except for 
the following differences. OSA performs detection rarely, and the PU is stable for the duration of OSA 
network operation, i.e. it is either transmitting constantly on a channel or stays idle. In other words quiet 
period occurs for multiple time slots, see Fig. |l(b)[ Also, since the PU is considered stable on every 
channel we do not consider all types of OSA MAC protocols introduced for the microscopic model. 
Instead we use classical DCC and HCC models proposed in f34| with the corrections of |fT5l accounting 
for the incomplete transition probability calculations whenever OSA network occupied all PU channels 
and new connection was established on the control channel. 

B. Formal Problem Description 

To compute the maximum throughput for different combinations of protocols and models, we define an 
optimization problem. The objective is the OSA network link layer throughput Rt. Therefore, considering 
the regulatory constraints given above we need to 

maximize Rt = £,R subject to = Pd,mm: td < id,max, (1) 

where td is the detection time, i.e. the time to process whole detection operation as described in 
Section IIII-D[ R is the steady state link layer throughput without sensing and switching overhead, which 
will be computed in Section ITVl and 

tt-tg-tp microscopic model and DCC with channel switching, 



td.„ax-t, 



macroscopic model, (2) 



^V^i otherwise, 

It 



10 



TABLE I 

Summary of Variables and Functions Used in the Paper Grouped by Relation 



Parameter 


Description 


Unit 


C, b 


channel througiiput, and channel bandwidth 


(M)bps, (M)Hz 


M, (Md), a, N 


total (effective) number of channels, fraction of M, and number of users 




Rt, R 


throughput: total, excluding sensing overhead 


(M)bps 


Qp Pc 


PU activity level, and PU activity seen by SU 




1/9, d 


average packet size 


slot, bit 


P 


SU channel access probability 


_ 


^t) iu^ iqri ir 


length: total slot, slot data, switching time, quiet time, sensing, reporting time 


(^)s 


te-, ta 


time to sense, and transmit sensing information by one node 


(m)s 




detection delay, and maximum detection delay 


(p)s 




slot overhead ratio 




e, K, 6, 7 


time-bandwidth product, cooperation level, detection threshold, PU signal level 


(M)Hzs, — , dB, dB 


r(-),r(-,-) 


complete and incomplete gamma function 







function supporting computation of optimal 9 


_ 


Pd, Pd,i, pii, pii 


detection probability: total, group, individual, individual with error 




Pf' PS,i' PW, PW 


false alarm probability: 




Px, Px 


supporting function to compute probability of detection and false alarm 






minimum probability of detection 




Pe 


probability of error during transmission 






number of: channels sensed per group, users in group, groups 




mi,i, m2,i, Vi 


supporting expressions to derive fhr,i for TTDMA 


bit 


fhr, fhr,i 


network, and group average number of bits to report 


bit 


fhs 


number of sensing cycles 







number of zeros, and number of bits until the end of report for TTDMA 


bit 


t, i, j 


time instance, and dummy variables 




«S, Sm 


state space of Markov chain, and size of <S 




k, I, m 


number of: channels used by SU, used by PU, SU coimections (previous slot) 




X, y, z 


(current slot) 




Xu Yt, Zt 


(at t) 




"^xyzt Pkl{7n)\xy{z) 


steady state, and transition probability 




'^x^y-j 


functions supporting Pki\xy 




^x,y^ J^x^y 


total PU arrival probability for no buffering, and buffering case 






termination, and arrangement probability 




c(l) 


functions supporting 




N 


set of natural numbers 




I, U 


vector of active users, and free channels 




Ii,t, Ili,t, Uj,t 


node i index, priority of node i, and free channel j index at t 




a, b, c {d, e, /) 


users' (channel) enumerators 






indication function in channel switching 




l/qu 1/qt, (1/fe) 


SU packet length i, and at time ti (t2) 


slot 
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TABLE 11 

Summary of Abbreviations Used in the Paper Listed Alphabetically 



Abbreviation 


Explanation 


DCC 


Dedicated Control Channel 


HCC 


Hopping Control Channel 


MAC 


Medium Access Control 


OSA 


Opportunistic Spectrum Access 


PHY 


Physical Layer 


PU 


Primary User 


SPCC 


Split Phase Control Channel 


SSMA 


Single Slot Multiple Access 


SU 


Secondary User 


TDMA 


Time Division Multiple Access 


TTDMA 


Truncated TDMA 



where 1 — ^ is the sensing and switching overhead, see also Fig. [T] Note that in ([U is itself affected 
by pj, as it will be shown in Section |IV] Also note that tp is removed from second condition of Q since 
the switching time is negligible in comparison to inter-sensing time. 

III. Layered Model of Spectrum Sensing Analysis 

To design the spectrum sensing, we follow the approach of iTJ in which the spectrum sensing process 
is handled jointly by (i) the sensing radio, (ii) the sensing PHY, and (iii) the sensing MAC. Using this 
layered model we can compare existing approaches to spectrum sensing and choose the best sensing 
architecture in a systematic way. Since the parameters of the design framework in ([T]) are determined by 
the choices of individual layers, we describe and parametrize each layer of the spectrum sensing, later 
describing cross-layer parameters. 

A. Sensing Radio 

The sensing radio scans the PU spectrum and passes the spectrum sensing result to the sensing PHY 
for analysis. The sensing radio banwidth is given as aMb, where a is a ratio of the bandwidth of the 
sensing radio to the total PU bandwidth and 5 MHz is the bandwidth of each PU channej^. With a > 1/M 
node can sense multiple channels at once. However the cost of such wideband sensing radio increases. 

^Note that C used later in calculating MAC throughput is an average throughput obtained using a certain modulation over a 
channel with bandwidth b. 



12 



B. Sensing PHY 

The sensing PHY analyzes the measurements from the sensing radio to determine if a PU is present 
in a channel. Independent of the sensing algorithm, such as energy detection, matched filter detection 
or feature detection BHll . B9l . there exists a common set of parameters for the sensing PHY: (i) time 
to observe the channel by one node /"S, (ii) the PU signal to noise ratio detection threshold 6, and 
(iii) a transmit time of one bit of sensing information ta = 1/C/is. We denote conditional probability 
of sensing result pij,i,j G {0,1}, where j = I denotes PU presence and j = otherwise, and i = I 
indicates the detection result of PU being busy and i = otherwise. Observe that piQ = 1 — pqq and 
pii = 1 - poi. 

As noted in Section III-AI we consider energy detection as the PU detection algorithm since it does 
not require a priori information of the PU signal. For this detection method in Rayleigh plus Additive 
White Gaussian Noise channel pio is given as ifTSl Eq. (1)] 



Pio = V, ' , (3) 



r(e) 



and pn IBl Eq. (3)] 



Pii = e 2 




e-7 



(4) 



h=0 

where r(-) and r(-, •) are complete and incomplete Gamma functions, respectively, and e = [teaMb\ is 
a time-bandwidth product. By defining Ge{9) = Pm and 9 = G7^(pio), we can derive pn as a function 
of pio and tg. 

C. Sensing MAC 

The sensing MAC is a process responsible for sensing multiple channels, sharing sensing results with 
other users, and making a final decision on the PU presence. Because of the vast number of possibilities 
for sensing MAC algorithms it is hard to find a general set of parameters. Instead, we derive cross- 
layer parameters for a specific option of the sensing MAC. This methodology can be applied to any 
new sensing MAC scheme. We now introduce classifications which will be used in the derivation of 
cross-layer parameters. 

7 ) Sensing Strategy for Grouping Channels and Users: Each SU has to determine which channel 
should be sensed among the M channels. To reduce sensing and reporting overhead, OSA system can 
divide users and channels into Ug sub-groups jSOll . Sub-group i G {1, • • • ,ng} is formed by n„ j users 
who should sense rris j channels to make a final decision cooperatively. Assume that all users are equally 
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divided into groups then rus^i € {[M/n^J, [M/n^]} and nu,i G {[N/ng\, iN/ng]}. Note that for 
M/ug G N and N/ug G N all sub-groups have the same nu,i = N/ug and m^^j = M/ug for all i. Given 
N and M, if is small, more users are in a group and the collaboration gain increases, but at the same 
time more channels must be sensed, which results in more time overhead for sensing. For large Ug, this 
relation is opposite. 

2) Combining Scheme: By combining sensing results of other users, a OS A network makes a more 
reliable decision on PU state. As considered in |[T3l . lISTI . we will take k as a design parameter for the 
sensing MAC and find an optimum value to maximize the performance. Note that for the case of user 
cooperation if k = 1, the combining logic becomes the "or" rule |fT9l Sec. 3.2], [42, Sec. III-C] and if 
K = N, it becomes the "and" rule. 

3) Multiple Access for Measurement Reporting: To transmit sensing results of multiple users through 
the shared media, a multiple access scheme is needed. Note that this multiple access scheme is only 
for the reporting process, different from the multiple access for data transfer. We consider the following 
approaches. 

a) Time Division Multiple Access (TDMA): This is a static and well-organized multiple access 
scheme for which a designated one bit slot for sensing report transmission is assigned to each user |43 |, 

m. 

b) TTDMA: In TDMA, when the SU receives all the reporting bits from other users the SU makes 
a final decision of presence of PU on the channel. However, in OSA network using TTDMA SUs may 
not need to wait until receiving the last reporting bit, because for the "k out of A^" rule, a reporting 
operation can stop as soon as k one bits denoting PU presence are received. This sensing MAC aims at 
reducing the reporting overhead, but unfortunately we have not seen any paper proposing and discussing 
TTDMA. 

c) Single Slot Multiple Access (SSMA): For this scheme, known also as the boosting protocol |[52l . 
only one bit slot is assigned for reporting and all SUs use this slot as a common reporting period. Any SU 
that detects a PU transmits one bit in the common designated slot. Otherwise, a user does not transmit any 
bit in the designated slot. Then, reporting bits from SUs who detect a PU are overlapped and as a result 
all power of the slot is summed up. By measuring the power in the designated slot, a SU can determine 
whether the primary user exists or not. We assume perfect power control and perfect synchronization. 
Even though this may not be practical, because carrier frequency or the phase offset cannot be avoided 
in real systems, this scheme serves as an upper bound for sensing MAC performance. For the analysis 
of SSMA in isolation but in a more realistic physical layer conditions the reader is referred to |[53l . |[54l . 
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D. Cross-Layer Parameters 

Considering the combined impact of the individual layers, we derive cross-layer parameters in the 
framework as described in ([T|). More specifically these are tq and t^, derived as a function of individual 
parameters and pf, and p^, denoting final network-wide probabilities of false alarm and detection, 
respectively. 

1) Detection Time td and Quiet Time tq.- Detection time td is defined as the time duration from the 
point that a SU starts to sense, to the point that a SU makes a final decision on PU presence. Regardless 
of the data transfer and spectrum sensing time overlap, the final detection decision is made only after 
combining the sensing group's reported information |[55l . Thus td is the time from the start of the sensing 
phase to the end of the reporting phase, i.e. td = ts + tj.. 

Since the data transfer may not be possible during sensing or reporting phases tq < td, depending on 
the approach. When spectrum sensing and data transfer are divided in time division manner tq = tg + tj.. 
Note that three other methods sharing the same problem are possible (they will not be considered in the 
remainder of the paper): (i) simultaneous reporting and data, which can be implemented by using the 
separate channel as in [56|, for which tq = ts, (ii) simultaneous sensing and data, implemented by using 
the frequency hopping method as in ||57l . for which tq = tj., and (iii) simultaneous sensing, reporting, 
and data for which tq = 0. Conceptually, simultaneous sensing, reporting, and data transfer is possible 
and seems most efficient but we have not found any implementation of it in the literature. Note that 
in order to implement simultaneous sensing and transmission at least two radio front ends are needed, 
which increases the total cost of the device. 

Define ifig as the number of individual sensing events to complete sensing operation and ifij. as the 
average number of bits to report. Then the sensing time and the reporting time can be calculated as 
ts = fhste and tr = mrta- Note that fhs is affected by the bandwidth of the sensing radio because it can 
scan multiple channels at once if the bandwidth of the sensing radio is wide. For the case that the sensing 
radio is narrower than the bandwidth to sense, i.e. a < maxjms^i, • • • , ms^Ug}/]^, we assume that a SU 
monitors all channels by sequential sensing |[33l . because the reporting phase should be synchronized after 
all SUs finish the sensing phase. With this assumption tjIs — |^niax{?7is i, • • • ^Tfis^n, 

^ } / qM] , because 

even though the bandwidth to sense is less than that of the sensing radio it still needs one sensing cycle 
to get information. For fhj., because there are Ug groups in a OS A system, fhj. = Y^^=i "^r.i where j 
depends on the multiple access schemes for reporting, which we compute below. 

a) TDMA: All nu,i users should transmit the sensing results of nis^i channels. Thus, fhr^i = Uu^ims^i- 
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b) TTDMA: For k < j/2, if k of ones are received, the reporting process will end. We introduce 
a variable 6 which is the number of bits when the reporting process finishes. Thus there should be k — 1 
of ones within 6 — 1 bits and then 5-th bit should be one. Because the range of 6 is from k to n„ j, the 
average number of bits for this condition is derived as 

"^1.^ = E (! 1 1) {(1 - Mo^Pw + QpSptrPn} ■ (5) 

Moreover, if the number of received zeros, denoting PU absence, equals to n„ j — k + 1, the reporting 
process will stop because even if the remaining bits are all one, the number of ones must be less than 
K. Then the reporting process stops at 5-th bit if 5 — riu^i + k — 1 bits of one are received within 5 — 1 
bits and zero is received at 5-th bit. The range of 5 is from j — k + 1 to n„ j, and thus the average 
number of bits for this condition is 

"^2,* = U _ M {(1 - QpWqoPw''' + (IphmPil'''] , (6) 

5=Vi 

where Vi = n^^i — k + 1. Therefore because there are m^^j channels to sense in a group i, fhr^i = 
ms,i{mi^i + m2,j). 

For the case k > n„ j/2, m-i j is calculated by counting zeros and m2^i by counting ones. Thus, we 
use fhr^i = ms^i{mi^i + m2,j) again, by replacing k with n„^j — k + 1, poo with piQ and poi with pu. 

Because we assumed so far that k is known to each node in the network, OSA nodes know when to 
stop reporting measurements and start data communication without being instructed by external parties. 
For comparison we analyze another type of TTDMA, denoted as kTTDMA, where a cluster head node 
makes a decision to stop the reporting phase in the OSA network. For example, this approach may be 
necessary if the k value is updated in real time. In the worst case scenario this approach requires two 
bits to be reported by the SU, i.e. one for sending sensing data and one for an acknowledgment from 
the cluster head to report. Then ^ and ^ need to be modified by multiplying them by two. 

c) SSMA: For this scheme, we need only one bit per channel for reporting. Thus j = j. 

2) Total False Alarm Probability pf and Detection Probability p^: Final probabilities pf and pd are 
obtained after cooperation, and thus affected by the sensing MAC and sensing PHY. Because each sub- 
group has a different number of users and channels to sense, we have 

^ rig 

Pf = —y^Pf,u (7) 
" 1=1 

where j is the probability of false alarm of sub-group i. Using (O we can also derive p^i by substituting 
Pf^i with pd^i, i.e. probability of detection of sub-group i. The definitions of and for each protocol 
are as follows. 
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a ) TDMA: For this protocol p f^i is derived as 



5=K ^ 



''U,i \ i—S 

g I PioPoo ' 



(8) 



where px = (1 — Pe)Px +Pe(l — Px) for Px G {piO;Poo}> while p^i is derived from ([8]l by substituting 
Pio with pii and poo with poi- 

TTDMA: In this case SU does not need to receive j bits to make a final decision because 
the reporting phase is ended when the number of ones is n. To derive j for this case, we introduce a 
variable /3 denoting the number of zeros. Then total number of reporting bits is k + /3 if the last bit is 
one because otherwise reporting phase will end at less than k + /3 bits. Therefore, there should be /3 of 
zeros in K + /3 — 1 bits and K-th bit should be one. Because (3 can vary from to riu^i — k 



Finally j is obtained from ^ by substituting piQ with pn and poo with poi. 

c) SSMA: Obviously, the process of the reporting information for SSMA is the same as for TDMA. 
Therefore pf^i and p^ j are defined the same as for TDMA. 

IV. Multichannel OSA MAC Protocol Analysis 

In this section we present the analysis of throughput R for all considered combinations of MAC 
protocol architectures. As noted in Section II-CI we propose a set of new multichannel MAC protocols 
for OSA. We will first describe their operation, later presenting the analysis framework. 

A. Description of New Multichannel MAC protocols for OSA 

We consider two major groups of MAC protocols for OSA: (i) those enabling buffering of the SU 
connections preempted by the PU arrival, and (ii) those enabling switching of the SU connections to a 
vacant channel when preempted. In the former group, when the PU arrives the existing SU connection will 
pause at the time of preemption and resume on the same channel as soon as the PU goes idle. We assume 
that the SU always waits for the PU to finish its transmission. The case where the buffered SU connection 
expires after a predefined time, not analyzed here, is presented in ll22l for the centralized network. We 
do not consider any channel reservation schemes for potential SU connections to be buffered |[25l . When 
buffering is not possible, the preempted SU connection is considered as lost and a new connection 
must be established on the control channel. In the latter group, when the PU arrives the existing SU 
connection will look for a new empty channel, to continue transmission. If such a channel cannot be 




(9) 
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found the connection is lost. Without channel switching, the exiting SU connection is lost as soon as the 
PU preempts the channel. 

Obviously we can have four combinations of these groups for OSA MAC, which have all been 
considered in the analysis: (i) with no buffering and no channel switching |30] scheme denoted as BqSo, 
where SU connections preempted by PU are lost; (ii) with no buffering and channel switching |[24l . |[25l . 
||26l denoted as BqSi, where SU connections preempted by PU switch to a free channel and connections 
that cannot find a free channel are blocked; (iii) with buffering and no channel switching ifTSl . ll22l . 
||23l denoted as BiSq, where SU connections preempted by PU are being suspended from the moment 
of preemption until PU releases the channel; and (iv) with buffering and channel switching 1581 denoted 
as BiSi, where SU connections preempted by the PU first looks for free channels, and when no free 
channels are found, the connections are being buffered until PU leaves the channel. The detailed procedure 
of distributed channel selection on the event of switching will be described in Section IIV-C2I Recall that 
the works referred above consider an OSA network with centralized channel management, providing only 
the upper bound on OSA network performance. 

B. Multichannel MAC for OSA Analysis: Preliminaries 

Usually, to compute the throughput of most non-OSA network it is assumed that a Markov chain 
characterizes the network state defined as the current number of connections used for data transfer |[T5l 
Sec. III-C], ||34l Sec. 3]. The state transition probability depends only on the network users' traffic 
characteristics. However, in the OSA system the SU data transfer connections can be terminated or 
delayed if a PU is detected on their channel, and thus the traffic generated by PU also affects the 
state transition. Moreover, with connection buffering enabled, when a PU is detected on the channel the 
SU does not terminate its connection but rather waits until the PU goes idle. Thus the OSA network 
throughput is influenced by the number of channels that are actually utilized by the SUs rather than 
solely by the number of SU data transfer connections. 

We propose a three dimensional Markov chain of which the state vector is given as (Xt, Yt,Zt), where 
Xt, Zt G 5 = {0, 1, . . . , min([iV/2j , Md)] and G {0, 1, . . . , Md], where M^, = M - 1 for DCC and 
Mn = M for HCC. The elements of the state vector are: (i) Xt denoting the number of channels that 
are actually utilized by the SUs at time t, (ii) Yt denoting the number of channels on which the PU is 
detected at time t, and (iii) Zt denoting the number of connections for the data transmission between the 
OSA users at time t. This distinction allows to compute the exact channel utilization, contrary to |[22l 
Sec. Ill] where buffered SU connections were also considered to be utilizing the PU channels. 
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Considering a real OSA system, there are conditions that quaUfy vaUd states. With SU connection 
buffering-enabled MAC protocols for OSA, the number of connections cannot be less than the number 
of channels utilized by SUs, i.e. Xt < Zf. Additionally, SUs do not pause transmissions over unoccupied 
channels. Therefore, the number of SU connections not utilizing a channel cannot exceed the number of 
channels occupied by PUs, i.e. Zt — Xt < Yt or Zt < Xt + Yf. Finally, the sum of the channels utilized 
by PUs and the SUs cannot be greater than Md, i.e. Xt + Yt < Mo- By combining these conditions we 
can compactly write them as 

0< Xt<Zt< Xt + Yt <Md. (10) 

When connection buffering is disabled the number of SU connections must be the same as the number 
of channels utilized by SUs, i.e. Xt = Zt. Therefore, for non-buffering SU connection OSA MAC 
protocols {Xt,Yt,Zt = Xt) {Xt,Yt). 

For the microscopic case the average channel throughput, excluding switching and sensing overhead, 
is computed as 

R = C^^^X7r,y,, (11) 

x=0 y=0 z=0 

where Sm = max{5} and the steady-state probability Hxyz is given by 

TT^y:, = lim Fr{Xt = x,Yt = y, Zt = z), (12) 

and the state transition probabilities to compute ([T2l ) will be derived in the subsequent section, uniquely 
for each OSA multichannel MAC protocol. 

Finally, for the macroscopic case the average channel throughput, excluding switching and sensing 
overhead, is computed as 

R = {qp{l - pd) + (1 - qp){l - Pf)}RcC, (13) 

where Rc = X]i=i '^'^i ^^'^ '^i ^ solution to a steady state Markov chain given by ifTSl Eq. (13)]. 
Since the macroscopic model assumes no PU activity in each time slot, SU connection buffering and 
switching is not needed. Note that contrary to the incorrect assumptions of |[T5l Eq. (12)], ll34l Eq. (7) 
and (9)] we compute R in ([TT] ) and (fT3l ) taking all the channels into account, irrespective of the type 
of OSA MAC. This is because models of ifTSl . |[34l considered only data channels for the throughput 
investigation in DCC in the final calculation stage, assuming that no data traffic is being transmitted on 
control channel. However, the utilization must be computed over all channels, irrespective of whether 
one channel transmitted only control data or not. 
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C. Derivation of State Transition Probabilities for the Microscopic Model 
We denote the state transition probability as 

Pxyz\kim = Pr{Xt = x,Yt = y, Zt = z\Xt-i = k, Yt-i = I, Zt-i 



m) 



(14) 



Note that changes in Xt and Zt depend on the detection of the PU. In addition, changes in Zt depend on 
OS A traffic characteristics such as the packet generation probability p and the average packet length 1/q. 
Also, note that the steady state probability vector vr containing all possible steady state probabilities iTxyz 
is derived by solving vr = vrP, where entries of right stochastic matrix P are defined as (fT4l ) knowing 



x,y,z "xyz 



1. 



As a parameter to model PU state, pc denotes the probability that a OSA network collectively detects 



a PU channel as occupiecj^, i.e. 



Pc 



QpPd + (1 - qp)pf- 



(15) 



.(i) 



.(i) 



We introduce two supporting functions. First, we denote Tjf as the probability of termination of j 
SU connections at time t given that k channels are utilized by the OSA network at time t — 1, which is 
derived as (H Eq. (2)] 

(16) 

0, otherwise. 

V 

Note that k in denotes the number of channels utilized by OSA network rather than the number of 
SU connections because only active connections can be terminated at the next time slot. And second, we 
denote Sm as the probability of j SU successful new connections at time t, given m connections were 
active at time t — 1. We need to modify the definition of 5"^^ given in ll34l Eq. (5) and (8)] considering 
PU detection on the control channel. If a PU is detected on a control channel, an SU connection cannot 
be generated because there is no chance to acquire a data channel. We then have 1,15. Eq. (17)] 



J 



r.(l) N-2m-l Mp-m 
•-"m AT-l M ' J 

1 - 9^') 



0, 



1 (DCC), 
1 (HCC), 

i = o, 

otherwise, 



(17) 



^Note that, contrary to 1151 Eq. (8)], we do not consider packet capture effects in the definition of pc as the packet capture 
model proposed in 1151 was an approximation. 
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where 



Sm\ PU free control channel, DCC only, 

(18) 

(1 — Pc) Sm ^ , Otherwise , 



and = {N — 2m)p{l—p)^ ^™ ^. Again, note that for the SU connection buffering protocols the sub- 
parameter m of Sm is not the number of channels utilized by SUs, but the number of SU connections. 
This is because we assume that a SU that has a connection but pauses data transmission due to the 
PU presence does not try to make another connection. We can now derive the transition probabilities 
individually for all four different OSA MAC protocols. 

1) Case BqSq: Recall that for non-buffering OSA MAC protocols Zt = Xt. Thus Pki\xy is defined as 
([141 ) without Zt- Because it is assumed that no more than one connection can be generated in one time 
slot, it is impossible to transit from k connections at time t — ltox>A; + l connections at time t. The 
state transition probability for this condition is 0. 

For X = k + 1 only one SU connection is created, no current connection is terminated and y PUs can 
appear on the channels that are not utilized by SU. Thus, from Mr, — x channels, a PU appears on y 
channels, so cases of PU appearances are possible. 

Now, consider the case x < k + 1. When a SU data connection is terminated, there can be two possible 
reasons: (i) a SU completes its transmission, or (ii) a PU is detected on a channel that is assigned to a 
SU for data transmission before sensing. The former was analyzed |[34', Sec. 3]. To model the latter, we 
introduce variable i denoting the number of channels that were reserved for SU data transmission before 
sensing but cannot be utilized due to PU detection. We have the following observation. 

Observation 1: For multichannel OSA MAC without SU connection buffering or channel switching 
the number of PU appearance combinations is (^^*) ■ 

Proof: When the OSA network detects PU on i channels from x + i channels that are going to be 
utilized by SUs before sensing, there can be (^^*) possible combinations for the PU appearance on the 
channels. For the remaining Mr, — x — i channels, y — i channels should be occupied by the PU because 
the total number of channels in which a PU is detected at time t should be y. Thus there are (^^"Z^ 
possible combinations for PU appearance on unassigned channels. □ 

For the SU traffic generation in the case of x < A; + 1, there are two possible cases: (i) no connection 
is created and k — x + i connections are terminated, and (ii) one connection is created and k — x — i + 1 
connections are terminated. Recall that k connections at t — 1 are changed to x + i connections at t 
before sensing. Also note that i G [0, 1, . . . , min (sm — y)\, where Sm — x is the number of possible 
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SU connections that can be terminated by PU appearance. By summing over all possible i we can compute 
the transition probability for the case x < A; + 1. 

In addition, we need to discuss the edge statcl which considers two cases: (i) no more channels are 
available, either utilized by SUs or PUs, and (ii) all possible SU connections are establishecjfl which we 
denote as "full connection state". For the transition from full connection state to edge state, we have to 
consider the case that one new connection is generated while any existing connection is not terminated, 
which means a trial for the new connection by the free SU is not established because there already exists 
all possible connections. 

Writing all conditions compactly, denote the indicator for the edge state 



<-x,y 



1, x + y = Md or X = Sm, 
0, otherwise, 



(19) 



and define Pxl 



i^V) C'^^y-l *)?'c(l - PcY^°^y, the complete state transition probability is given as 



Pxy\kl 



0, 



7^(0) c^(l) p(0) 



X > k+1 
X = k + 1, 

X < k + 1, k < Sm or Ix ' 



i=0 



(20) 



E/ ^{k-X-i) ^(0) _^ rj^ik-X-i + l) p{i) 

i=0 



x,y 



I 7^(0) p(0) 
^ ^k '^k ^0,y ' 

where im = niin(sm — x, y). 

2) Case BqSi: Although in the SU connection non-switching case both DCC and HCC can be 
considered, only DCC will be able to perform switching without any additional control data exchange, 
which we prove formally. 

Before going into detail of the derivation note that for the class of OSA MAC protocols with a dedicated 
control channel every node can follow the connection arrangement of the entire network. Because the 
dedicated control channel is continuously monitored by all network nodes via a separate front end ll34l 
Sec. 2.2], each node can learn the overall network configuration. Note that this also applies to Split Phase 



''As shown in [15. Sec. III-D] the edge state was not considered in |34 Eq. (6)], which resulted in an incorrect model. 

^If Sm = Md there can be many free SUs but no channel available. On the other hand, if Sm ~ \_N/2\ there can be one 



free SU for even A'^ or no free SU for odd N. 
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Control Channel MAC (SPCC) IMl Sec. 2.4], HI Fig. 2(c)] as well, since SPCC has a dedicated control 
channel phase. For HCC, as well as Multiple Rendezvous Control Channel |[34l it is impossible for a 
single node to learn the whole network connection arrangement since each sender receiver pair cannot 
listen to others while following its own hopping sequence. We now present the following proof. 

Theorem 1: Channel switching in DCC can be performed without any additional control message 
exchange. 

Proof: We prove this by showing a possible distributed channel switching process. Following earlier 
observation, in DCC each node can trace the connection arrangement of others, i.e. which channel has 
been reserved by a sender receiver pair. To distribute the switching events equally among SUs each SU 
computes the priority level as 

ni,t = iii^t-i + ip, (21) 

where 

1, preemption by PU, 

(22) 

0, otherwise, 



and Ilj ( is the priority level of SU i at time t. For Ilj o ^ N the priority is a MAC address of the SU, 
transformed into a real number for each SU by a network-wide known function. Now, having a set of 
priorities of all communicating node pairs, the OSA network is able to select, in a distributed manner, a 
new set of communication channels upon PU arrival as 

{la^uh,U Ic,t} ^ {Ud,t, Ue,t, t//,t}, (23) 

V ' V ' 

I u 

where |I| = |U| = Md — Xt — Yt, ^ is the mapping operator denoting process of switching active SU 
connection i to free channel j, li^t denotes index of communicating SUs (transmitters) at time t, where 
Iia,t > ^b,t > • • • > and Uj^t denotes free channel with index j at t. □ 
Note that existing connections that have not been mapped to a channel are considered blocked. Also note 
that algorithm given in Theorem [T] connections are preempted randomly with equal probability by PU. 
Since new SU connections are also assumed to use new channels randomly with equal probability, each 
SU connection is blocked with uniform probability. 

To enable SU connection switching in HCC one way is to augment it with a separate radio front end 
which would follow the hopping sequences and control data exchange of the OSA network. Obviously 
this increases the cost of hardware and contradicts the idea of HCC, where all channels should be used 
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for data communication. Therefore while evaluating OSA MAC protocols in Section IV-BI we will not 
consider SU connection switching for HCC. 

We now define the state transition probability p^^i^; for the considered OSA MAC protocol. Because 
X > /c + 1 is infeasible, the state transition probability for x > A; + 1 equals to zero. For x = k + 1, y 
PUs can appear on any of Mu channels because even though a PU is detected, the SUs can still transmit 
data by switching to the idle channels and the possible number of PU appearances is (^f°). Note that the 



possible number of PU appearances in the case BqSi is always 



fMr 



), even for the edge state, because 



the data channel can be changed by switching to a vacant channel after the PU detection. Because it is 
impossible to create more than one new connection at a time, the OSA connection creation probabilities 
for x = + 1 are the same as in (l20l) . i.e. xjf^ 

For X < + 1, if SUs are not in an full connection state, there are two cases for the OSA traffic 
generation: (i) no connection is created and k — x connections are terminated, and (ii) one connection 
is created and k — x + I connections are terminated. On the other hand, for the state transition to the 
edge state we use variable i, just like in the case BqSo, to derive the probabilities because connections 
can be terminated by a PU in the full connection state state. Furthermore, for the transition from the full 
connection state to the full connection state, we should take into account Tjf^ S^'^ again. 

With all conditions the state transition probabilities are denoted compactly as 



Pxy\kl 



rpiO) p(0) 

-^k '^k ^o,y 

(x-k) ^,(0) 



(x-k+1) p(0) 



n{k-X~i) n{0) rp{k-X-i+l) c-(l) 
■^fc ^-'fe '^k 



Pr 



(0) 



0,1/' 



X > k + 1, 

X = k + 1, 

X < k + 1, 1^ y = 0, 

X <. k -\- 1, k <^ S-mt '^x,y 



1, 



(24) 



i=0 



i=0 



-,{k-x-i) n{0) 



(k-x-i+l) c.(l) 



(0) 



0,y 



1. 



, 7^(0) c;.(l) p(0) 
"•"^A: '^k ^0,y> 

3) Case BiSq: Before we discuss this case we present the following observation, which implicates the 
design of simulation models and derivation of Pxyz\kim for SU connection buffering MAC protocols. 

Observation 2: For all SU connection buffering OSA MAC protocols the same average link level 
throughput results from creating a brand new connection or resuming a previously preempted and buffered 
connection on the arrival of PU on a channel. 
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Proof: Due to the memoryless property of the geometric distribution 

Pr(l/(7i > l/qt, + l/qu\l/qi > 1/fc) = Pr(l/% > l^J, (25) 

where 1/qi is the duration of connection i, l/qt^ is the connection length until time ti when it has been 
preempted by PU, and l/qt2 is the remaining length of the connection after SU resumes connection 
at time ^2- Since either a newly generated SU connection after resumption, or the remaining part of a 
preempted connection needs a new connection arrangement on the control channel, the number of slots 
occupied by each connection type is the same. □ 

Having Observation [2] we can derive transition probabilities. Because packet generation is affected by 
the number of connections, we use Zt to classify conditions to derive the state transition probabilities. 
Due to the assumption of a maximum number of one connection generation in one time slot, the state 
transition probability of the case of z > m + 1 is zero. 

For z < m + 1, a data connection is terminated only if the transmitting node completes its transmission 
without PU interruption. When a PU is detected in the channel, the SU temporarily pauses communication 
without terminating the connection. Among z SU connections, x connections actually utiUze channels 
for data transmission and z — x connections are paused due to PU detection. Thus (^^^) combinations 
of PU appearance are possible. At the same time, because the total number of channels occupied by PUs 
is y, the remaining y — z + x PUs should appear on — z idle channels. Thus there can be {y^_"^^l^ 
combinations for PU appearance on idle channels. 

The SU connection generation probability for z = m + 1 is Tjf^ Sm\ just like in the case BqSo. For 
z < m + 1 (i) no connection is generated and m — z connections are terminated, and (ii) one connection 
is generated and m — z + I connections are terminated. For z = m + 1, one connection is generated 
while no connections are terminated. 

For the transition between full connection states, i.e. the state transition from m = Sm io z = Sm, we 
should take into account the case that one connection is generated and no connections are terminated 
because there is no available resources for a new connection. 
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Finally, considering all these cases and defining Rx^y 
transition probability is given as 



0, 

rpifi) Al) „iz) 



Jiy-zZ)pUl-Pcr-y the state 



Pxyz\klm 



rp{m-z) q{0) rp(m—Z+l) 



R 



x,y, 



7^(0) q{0) , rpil) , ^(0) „{Z) 



'X,y, 



Z > 771 + 1, 

z = m + 1, 

z < m + l,m < Sm ot Z < Sm, 

Z = 171 = Sm.. 



(26) 



Note that this OS A MAC has been previously analyzed in lITSl . As it has been pointed out, the model 
proposed did not work well for the full range of parameters. This is due to the following. A Markov 
model has been derived for {Xt,Yt} (using unmodified transition probabilities of ll34l Eq. 6] used to 
calculate average throughput of networks based on non-OSA multichannel MAC protocols). With this 
limitation termination, the probability in |[T5l Eq. (14)], analogue to ( fT6l ). included an aggregated stream 
of PU and SU traffic, where PU traffic qp was later substracted from steady state channel utilization 
in |[T5l Eq. (10)], analogue to ([TT] ). The approximation of lITSl . although Markovian, worked reasonably 
well only for a moderate values of PU activity qp. 

4) Case BiSi: This OSA MAC from analysis perspective is the same as the buffering OSA MAC with 
no channel switching, except for the following two differences. First, if there is at least one idle channel, 
an SU that has a connection but does not utilize a channel cannot exist because this SU can switch to the 
idle channel. Formally, the state transition to the state of z / x and x -\-y < Mn and the state transition 
from the state m ^ k and k + I < Md are not possible. Second, in contrast to the non-switching OSA 
MAC, y PUs can appear in any of Mo channels because the SUs can switch to the idle channels in this 
option, the same as for BqSi case. Thus the possible number of cases of PU appearance is just (*^°). 
Therefore, replacing R^^^y with R'^y and adding conditions z ^ x, x + y < Md, m ^ k and k + l < Md 
to the condition z > m + 1 in ( [261 ) results in a complete definition of Pxyz\kim- For consistency we present 
this state transition probability as 
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5) Impact of Channel Error on the Throughput Calculations: All previous analysis were done under 
the assumption of the error-free channel. In this section we will briefly discuss the impact of channel 
error on the throughput calculations. 

Channel error impacts the throughput in two ways. First, error affects throughput when SU involved 
in a connection setup fails to receive a control message from the transmitter. As a result no connection is 
established. Second, error affects throughput when SU not associated with the current connection setup 
(which does not overhear the connection setup from others) collide later with other users, believing 
incorrectly it selected a free data channel for communication. The throughput of DCC is impacted by 
both effects. On the other hand, HCC is influenced only by the first effect. It is because HCC MAC 
protocol implementations do not posses a separate control channel. Thus no overhearing of connection 
setup is possible and users with HCC MAC protocol select data channel for communication automatically 
according to a pre-defined hopping sequence. Because of the prohibitive complexity of the analysis of 
the second effect, we focus on the first error case and the HCC. 

For HCC, the control channel is selected as one of the data channels by a hopping method. Thus, if we 
assume an error on the control channel, it is reasonable to consider the error on the data channel as well. 
For the control channel, if an error occurs, a connection fails to be established. Thus it is modeled by 
multiplying Sm by 1—pe, where p^ is a probability of error in the current time slot. For the data channel, 
different error handling strategies can be considered. We focus on the two following situations: i) case 
El denoting packet punctured by unrecovered errors and ii) case E2 denoting transmission termination 
on error. 

a) Case Ei: It can be assumed that when an error occurs on a time slot, the SU simply discards 
that time slot and resumes transmitting the remaining packet fragment from the next correct time slot. 
This is modeled by replacing the capacity C with C(l — pe)- 

b) Case E2: It can also be assumed that the connection terminates when an error occurs. Thus the 
probability that the packet finishes transmitting, q, should be replaced by g + (1 — q)pe- In addition, if 
the control channel hops to a channel which is being utilized for data transmission but error occurs, a 
new connection cannot be established. This is modeled by multiplying Sm by (1 — pe)^- 

V. Numerical Results 

We now present numerical results for our model. First, we present results independently for spectrum 
sensing and OSA MAC performance, in Section IV-AI and Section IV-BI respectively, for the microscopic 
case. Then in Section IV-CI we present the results of joint optimization of these two layers in the 
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microscopic and macroscopic context. Moreover, due to a vast combination of parameters to consider 
we have decided to follow the convention of lITSl . ll34l and focus on two general network setups (unless 
otherwise stated): (i) small scale network with M = 3, = 12, ci = 5 kB and (ii) large scale network 
with M = 12, N = 40, d = 20 kB. 

In this section we will also compare the analytical model of the sensing layer and OSA MAC protocols 
to simulation results. The simulations were developed with Matlab and reflect exactly the sensing models 
and MAC protocols presented in this paper. Simulation results for each system were obtained using the 
method of batch means for a 90% confidence interval. To evaluate the sensing protocols each batch 
contained 100 events and the whole simulation run was divided into 10 batches with no warm up phase. 
When simulating the OSA MAC protocols, each batch contained 1000 events while the whole simulation 
was divided into 100 batches with the warm up period equal of 100 events. 



A. Spectrum Sensing Architecture Performance 

For all possible combinations of sensing architectures we compute the probability of false alarm for 
a wide range of tq. For two networks considered we select a set of the following common parameters: 
tt = id,max = Ims, C = 1 Mbps, b = IMHz, qp = 0.1 (which approximately corresponds to the level 
of actual measured PU occupancy on the channel from 1591 Tab. 1]), 7 = — 5dB, and a = 1/M. In all 
sections, except for Section |V-A2[ we present results for an error-free channel. Note that for all results 
presented in this section simulation results for all protocols confirm the accuracy of the analytical model. 

1 ) Measurement Reporting Protocol Performance: The results are presented in Fig. |2] For the same 
Pd requirement, SSMA results in the lowest pf for each value of tq, while TDMA performs worst. The 
benefit of introducing TTDMA in comparison to TDMA is clearly visible for all network scenarios. 

The advantage of TTDMA and SSMA can be shown more clearly if we compare the results of different 
Pd = Pd,min requirements. We can observe that high detection requirement such as pd = 0.99 makes the 
performance worse, as generally known. However if TTDMA or SSMA is applied, the performance for 
Pd = 0.99 can be higher than that of TDMA for pd = 0.9. For example, in the range that tq < 50 /iS in 



Fig. [2(a)| SSMA for pd = 0.99 outperforms TDMA for pd = 0.9. Moreover, in Fig. [2(b)| for tq < 550 fis, 
SSMA and TTDMA for pd = 0.99 outperforms TDMA for pd = 0.9. 

It is important to note that kTTDMA performs worse than the rest of the protocols. It is due to 
excessive delay caused by instant acknowledgment of reporting result to the cluster head node. Note 
that kTTDMA is a lower bound for the operation of TTDMA. Also note that when TDMA needs to 
be equipped with acknowledgment function, as kTTDMA, its performance would be degraded the same 
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Fig. 2. Performance of different measurement reporting protocols as a function of tq and pd = Pd.min for (a) M = 3, N — 12, 
and (b) M — 12, N = 40. Common parameters: pe — 0, n = 1 ("or" rule), tt ~ fd,max = 1ms, C — 1Mbps, b = IMHz, 
qp = 0.1, 7 = -5dB, a = 1/M, and Ug = 1. 



way as TTDMA. Since we analyze static network with pre-set parameter values, e.g. k does not change 
over time, in the following sections we proceed with unmodified TTDMA only. 

2) Impact of Channel Errors during Reporting on PU Detection Performance: The results are presented 
in Fig. [3l For small and large scale network, and the same parameters as used in Section IV-All we have 
observed the probability of false alarm keeping detection probability pd constant for varying quiet time 
tq. First, it is obvious when comparing Fig. [2] (no channel error) and Fig. [3] (channel error) the impact 
of error is clearly visible, i.e. pf increases for every protocol. However, the relation between individual 
protocols is the same since error affects all protocols equally. Second, the effect of error on the small 
scale network is smaller than for the large scale network, compare Fig. 3(a) and Fig. |3(b)[ since the 
probability that SU will send a wrong report is larger for network with large number of nodes. Lastly, 
for small values of k probability of false alarm stabilizes and never reaches zero. However, large values 
of K reduce significantly the effect of channel errors. It is because with high k probability of making 
an error decreases rapidly. With 20% of nodes participating in the cooperative agreement on PU state, 
«; = 2 for small network and k = 8 for large scale network, effect of error is reduced almost to zero. 

3) Impact of Cooperation Level on PU Detection Performance: The results are presented in Fig. |4l 
We have selected TTDMA and set pd = pd,ram = 0.99 as a protocol for further investigation. We observe 
that for the small scale network, see Fig. |4(a)[ the performance for k = 2 is the best, while for the 
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Fig. 3. The effect of channel errors on the performance of different measurement reporting protocols as a function of tq and 
Pd = Pd.min = 0.99 for (a) M = 3, iV = 12, k = {1, 2}, and (b) M = 12, TV = 40, k = {1, 8}. All remaining pai'ameters ai'e 
the same as in Fig. |2] except for pe = 0.01 . 




large scale network, see Fig. |4(b)[ the best performance can be achieved when k = 8 or 16 if < 0.1. 
Based on this observation, we conclude that for given detection requirements, high detection rate of PU 
is obtained when k is well below the total number of SUs in the network. While for the considered setup 
optimal K w 20%, this value might be different for other network configurations . 
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4) Impact of Group Size on PU Detection Performance: The results are presented in Fig. |5] To contrast 
the impact of group size, we choose M = 12, N = 20 as the small scale network. We perform experiments 
only for the case when rhs equal for all groups, which means the number of groups is the divisor of M, 
i.e. Ug G {1,2,3,4,6,12}. 

An interesting observation is that the number of groups to achieve the best performance becomes larger 
as the number of users N increases. For the small scale network, see Fig. |5(a)[ the best performance is 
observed for = 2 or = 3, while for large scale network. Fig. |5(b) = 6 is the best. This is 



because for the large scale network, the reporting overhead caused by large number of users offsets the 
performance improvement achieved by large cooperation scale. 

5) Impact of k on PU Detection Performance: The results are presented in Fig. |6] For two network 
sizes, large and small, three sensing MAC protocols and fixed value of pd we vary k to see its impact on 
the sensing MAC protocol performance. We observe probability of false alarm as well as change in 
quite time tq. First we notice that varying k does not change probability of false alarm for any protocol, 
in both network configurations. Moreover the lowest probability of false alarm is obtained when small 
number of users agree on the PU state. The larger the channel number, the larger the range of n when 
network obtains the lowest probability of false alarm, compare Fig. |6(a)| and Fig. 6(b) The trends of 
for both network configurations are the same, since all protocols keep false alarm rate on the same level 
irrespective of the parameter change. In case of quiet time, TDM A and SSMA have qp constant and 
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Fig. 6. The effect of varying n on (a), (b) false alarm probability, and (c), (d) quiet time of different measurement reporting 
protocols as a function for pd = Pd.min and (a), (c) small scale networlc, M = N = 12, and (b), (d) large scale network, 
M = 12, N = 40. Rest of the parameters are the same as in Fig. |2] except for te = 20 ms. 



independent from At, which differs them from TTDMA whose operation strictly depends on the value of 
K considered. And again, when comparing Fig. |6(c)| and Fig. |6(d)| the optimal value of tq for TTDMA 
is in the same range as which proves the optimality of the design. 

B. OSA MAC Protocol Performance 

To evaluate the effectiveness of all proposed and analyzed MAC protocols we have fixed C = 1 Mbps, 
p = e^^/N, tq = tp = WO jjs, tt = td.min = 1 ms, Pd.min = Pd = 0.99, and pf = 0.1. Note that we 
do not relate p^ and pj with the actual spectrum sensing process at this moment (this will be done in 
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(a) (b) 

Fig. 7. Performance of OSA MAC protocols versus PU activity level for (a) i\f = 3, iV = 12, d = 5kB and (b) M = 12, 
A'' = 40, d = 20 kB. Common parameters: pe = 0, p = e^^/N, pd — pd.min ~ 0.99, pf — 0.1, tq = tp = 100 ^is, 
it = id, max ~ 1 ms, and C — 1 Mbps. 



Section IV-CI ). assuming that spectrum sensing layer is able to obtain such quality of detection. Again, as 
in Section IV-AI results are presented separately for error-free and error channel. 

1 ) Impact of PU Activity Level on OSA MAC Protocols: The results are presented in Fig. |7] We 
observe that PU activity degrades DCC and HCC for BqSo, irrespective of other network parameters. 
Their performances are comparable in this case. DCC and HCC performs best with B^Sq. The results 
show that the non-buffering OSA MAC protocols are very sensitive to qp where the greatest throughput 
decrease is visible at low ranges of PU activity. On the other hand, with connection buffering we observe 
a linear relation between qp and Rt. 

2) Impact of SU Packet Size on OSA MAC Protocols: The results are presented in Fig. [8] Obviously, 
for larger SU packet size, the OSA network is able to grab more capacity. However, when packets become 
excessively large the throughput saturates. It remains that with no buffering and no channel switching 
protocols obtain the lowest throughput, no matter what network setup is chosen. Interestingly, although 
intuitevely BiSi should obtain the highest channel utilization, it does not perform better than BiSq due 
to large switching time. With tp approaching zero, DCC BiSi would perform best, irrespective of the 
network setup as we discuss below. 

3 ) Impact of Switching Time on OSA MAC Protocols: The results are presented in Fig. |9l In this 
experiment, we verify that for small tp DCC BiSi outperforms DCC BiSq. However, there is no huge 
difference between their performances even at tp = 10 //s. This is because connection switching does not 
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Fig. 8. Performance of OSA MAC protocols versus packet size d for (a) M = 3, A*' = 12, and (b) M = 12, = 40. Rest 
of the parameters are the same as in Fig. [T] except for (jtp = 0.1. 
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Fig. 9. Performance of OSA MAC protocols versus channel switching time tp for (a) M — 3, N — 12, d = 5kB and (b) 
M = 12, A'' = 40, d = 20 kB. Rest of the parameters are the same as in Fig. [T] except for Qp = 0.1. 



seriously impact the data throughput for the network setups in which the number of channels is less than 
the number of possible connections, i.e. M < 2N. For this network setup, all channels are utilized in 
most of time, and therefore there may not exist many idle channels to switch. The performance of DCC 
BqSi is also improved for small tp, and we observe that DCC BqSi outperforms HCC BiSq for large 
scale network, see Fig. |9(b) 
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Fig. 10. Performance of OSA MAC protocols for ratio of number of SUs to number of PU channels for (a) M = 3, d = 5 kfi 
and (b) Af = 12, d = 20 kB. Rest of the parameters are the same as in Fig. [T] except for qp = 0.1. 



4) Relation Between Number of SUs and PU Channels: Finally we want to explore the relationship 
between the number of OSA network users and the number of available PU channels. The results are 
presented in Fig. [TOl 

With increasing ratio N/M we observe an increasing throughput, where at some point all protocols 
almost saturate. Again because of the high switching penalty DCC with BiSi is inferior to BiSq. For 
small scale network, as shown in Fig. |10(a)[ a separate comment is needed for HCC BiSq- For small 
N/M DCC with BiSi and BiSq obtains higher throughput than HCC BiSq. However, for high N/M 
HCC BiSo achieves the highest Rt of all protocols. For large scale network as shown in Fig. |10(b)[ 
comparing channel switching and buffering options we conclude that much more channel utilization is 
obtained by connection buffering than by channel switching alone when N/M > 1. 

Note that for all cases described in this section simulation results agrees with our analytical model. 
Comparing our model and analytical results of ifTSl for DCC BiSq, see Fig. |10(b)[ we observe that 
prior analysis overestimated the performance resulting in more than 2 Mbps difference at N/M = 1. 
Interestingly, if we consider the same set of parameters as in Section IV-B 1 1 then the model of |[T5l almost 
agrees with the model of our paper. Since the set of parameters that has been chosen in IV-B 11 are similar 
to fl5 \ we remark that the observations on the performance of this OSA MAC in ifTSl were reflecting 
the reality. 
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Fig. 11. Throughput of HCC OSA MAC as a function of packet size d for (a) M = 3, N = 12, pe = 0.1 and (b) M = 12, 
= 40, pe = 0.01 and two distinct error handling startegies. Rest of the parameters are the same as in Fig. [7] except for 
qp = 0.1. El and E2 denote error models described in Section llV-CSI Eo denotes the system with pe — 0. 



5) Impact of Channel Errors on the OSA Multichannel MAC Performance: To observe the impact of 
channel errors on the MAC protocol throughput we have set up the following experiment. For HCC and 
both network sizes, small and large, we have observed the average throughput for different SU packet 
lengths and channel error probabilities. The results are presented in Fig. [TT] For comparison in Fig. [TT] 
we present the system with no errors, denoted as Eq. We kept values of p^. realistic, not exceeding 1%. 
Obviously system with punctured errors Ei obtains much higher throughput than system E2, since more 
data can be potentially sent after one control packet exchange. Again, buffering allows to obtain higher 
throughput in comparison to non-buffered case, even with the data channel errors present. Note that 



system E2 is more prone to errors than Ei, observe Fig. |1 l(a)| and Fig. |1 1(b)] for BiSq Ei and BiSq E2. 

6) Impact of PU Channel Occupancy Distributions on the OSA Multichannel MAC Performance: All 
previous analyses were done under the assumption that traffic generated by SU and channel occupancy of 
PU can be described by the geometric process. This assumption holds generally either for SU traffic or 
PU channel occupancy statistics. For example, it has been shown recently in |[60l that geometric process 
constitutes more than 60% of the measured PU traffic in GSM 900 uplink, GSM 1800 downUnk, DECT 
and 2.4 GHz UNIT channels. It is important however to see the behavior of the considered data MAC 
protocols with other traffic distributions. Since the impact of different SU packet length distributions 
has been investigated in |[34l Sec. 5.2], concluding that comparable throughput of multichannel MAC 
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protocols is obtained, we focus on the impact of different PU traffic distributions on the OSA network 
performance. Due to vast number of combinations of protocol and traffic distributions we have narrowed 
our presentation to DCC and the following distributions: i) discrete uniform (denoted symbolically as 
U), ii) log-normal (denoted symbolically as L), and for comparison iii) geometric (denoted symbolically 
as E) used in the analysis. We have tested the protocol performance for different combinations of "on" 
and "off" times of PU activity. These were EE, LE, EL, LL (all possible combinations of "on" and 
"off" times obtained in ll60l Tab. 3 and Tab. 4]) and additionally EU, UU, where first and second letter 
denotes selected distribution for "on" and "off" times, respectively. Due to the complexity of the analysis 
we show only the simulation results using the same simulation method of batch means, with the same 
parameters as described at the beginning of Section |Vl 

The parameter of each distribution was selected such that the mean value of each distribution was 
equal to 1/pc for "on" time and 1 — 1/pc for "off" time. The uniform distribution has a non-continuous 
set of mean values, (at + an)/2, where ah, an S N denoting lower and upper limit of the distribution, 
respectively, which precludes existence of every mean on or off value for pc £ (Oi !)■ To solve that 
problem an continuous uniform distribution with required mean was used and rounded to the highest 
integer. This resulted in a slightly lower last peak in the probability mass function at a„ for 1/pc ^ N or 
1 — 1/pc ^ N. In case of log-normal distribution, because it is continuous, it was rounded it to the nearest 



where q = 1/pc, vi = — Pc)/p1 is the mean and variance of the resulting discretized log-normal 
distribution. Note that the variance of the used discretized log-normal distribution is equal to the variance 
of geometric distribution for the same mean value. The variance of resulting discretized uniform continuos 
distribution could not be equal to the variance of the geometric distribution due the reasons described 
earlier. 

The results are presented in Fig. [12] We focus on two network types, as indicated earlier: (i) large scale 
and (ii) small scale, with the assumed parameters as in Fig. |7] We select four values of qp for the clarity 
of the presentation. The most important observation is that irrespective of the considered distribution 
DCC obtains relatively the same throughput and the same relation between different protocol options 
exists as it was shown analytically in Fig. |7J If one wants to select the distribution combination with 
the highest throughput it would be LE and LL, while the throughput obtained being almost equal to the 
one obtained via analysis for the geometric distribution. The distribution with the lowest throughput is 
UU and EU, due to the difference of the second moment between the other two distributions for the on 
time. The difference in throughput between UU, EU and the remaining distributions is more visible for 



integer as well, with scale parameter a 




log f ^ + 1 j and location parameter /i = log 
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Fig. 12. Impact of different PU on and off times distributions on OSA DCC multiciiannel MAC performance; (a)-(d) small 
scale network, (e)-(h) large scale network, as described in as in Fig. [T] E, U, L denote geometric, uniform and log-normal 
distribution respectively, where the first and second parameter in the legend denotes off and on time, respectively. SN denotes 
small network and LN denotes large network. 



the large network. The most surprising result of this investigation is that any DCC MAC protocol option 
with buffering removes the impact of distribution type on the obtained performance, compare Fig. 1 12(a) 
and Fig. |12(c)[ or Fig. |12(e)| with Fig. |12(h)| for any value of qp. 

C. Performance of Joint Spectrum Sensing and OSA MAC Protocols 

Having results for spectrum sensing protocol and OSA MAC we join these two layers to form a 
complete OSA network stack. By means of exhaustive search we solve the optimization problem of ([Til. 
We will also investigate the set of parameters that maximize Rt for small and large scale network. 

We divide our analysis in macroscopic and microscopic case observing Rt for small scale network with 
M = 3, = 12, d = 5kB, and large scale network with M = 12, N = 40, d = 20 kB. For each case 
we select a set of spectrum sensing and OSA MAC protocols that are possible and, as we believe, most 
important to the research community. For a fixed set of parameters C = 1 Mbps, 6 = 1 MHz, p = /N, 
*(i,max = 1ms (microscopic case), td,ma.x = 2s (macroscopic case), a = l/M, tt = 1ms, Pd,mm = 0.99, 
7 = — 5dB, Qp = 0.1, and tp = 100 fis we leave k, t^, rig, and pf as optimization variables. 

1) Microscopic Model: Here we focus only on DCC protocol, since collaborative spectrum sensing 
is only possible via a PU free control channel, which is inefficient to accomplish with HCC. Also, 
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for sensing measurement dissemination we do not consider SSMA, which would be most difficult to 
implement in practice. The results are presented in Fig. \T3\ 



DCC BiSo with TTDMA is the best option, both for small scale and large scale network, see Fig. 13(a) 



and Fig. |13(b)l respectively. Because of relatively high switching time BiSi performs slightly worse than 
BiSo, for small and large scale network. DCC BqSo with TDMA is the worst protocol combination, which 
confirms earlier results from Section IV-AI and Section IV-BI Irrespective of network size it is always better 
to buffer SU connections preempted by PU than to look for vacant channels, compare again BiSq and 
BqSi in Fig. |13(a)| and Fig. |13(b)[ The difference between BqSo and BqSi is mostly visible for a large 
network scenario, see Fig. |13(b)[ since with a large number of channels there are more possibilities to 
look for empty channels. 

For all protocol combinations and both network sizes k = 2 maximizes throughput performance, see 
Fig. 13(a)[ Interestingly, network size dictates the size of a sensing group. For small scale network. 



= 1 is the optimal value, see Fig. |13(a)[ but for a large network Rt is maximized when Ug = 3 (for 
BqSo) and Ug = 4 (for the rest). We can conclude that with a small network it is better to involve all 
nodes in sensing, while for larger networks it is better to divide them into groups, which agrees with the 
observation from Section IV-A4I Moreover, we observe that the performance difference between TTDMA 
and TDMA is not as big as in Fig. |2] when parameters are optimized. 

The most interesting result is observed for pj. With the increase of protocol complexity false alarm 
increases as well. Also with an increase of pf, quiet time is decreasing. Because buffering and switching 
improves the performance, there can be more margin to design the spectrum sensing. 

2) Macroscopic Model: For the macroscopic model we explore both non-OS A DCC and HCC with 
TDMA and TTDMA as sensing protocols. The results are presented in Fig. [I4l 

DCC obtains higher throughput than HCC for a small scale network, and vice versa, compare Fig. |14(a)| 
and Fig. |14(b)[ respectively. This confirms the observations of I15i Fig. 3], 1341 Fig. 3]. Just Uke in 
Fig. 13(a)[ for small scale network k = 2 and Ug = 2 are the ones that maximize Rt. For the large 
scale network, however, k = 3 and = 3 is optimal for TDMA, and k = 4 and Ug = 4 for TTDMA. 
This means that for large networks it is beneficial to split the network into smaller groups. Again, this 
confirms our findings from Section IV-Cll For both network scenarios pj and is relatively the same 
for all protocols considered. 

Note that for the large scale network in the macroscopic model, an SU takes more time to detect 
a PU than in the microscopic model because large t^^max reduces the time overhead. The release of 
time restriction impacts the large scale network by requiring greater value of k to achieve the maximum 
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Fig. 13. Optimization result of the selected protocol combination with DCC for the microscopic model for (a) M = 3, N — 12, 

and (b) M = 12, N ^ 40. Common parameters: pe = 0, d = 5kB, C = 1Mbps, b = IMHz, p = e'~^/N, td,niB.x = 1ms, 
a = 1/M, tt = 1ms, pd.min = 0.99, 7 = -5dB, qp = 0.1, and tp = 100ms. 
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Fig. 14. Optimization result of the selected protocol combination for the macroscopic model for (a) M = 3, A'^ = 12, and (b) 
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throughput. 



VI. Conclusion 

We have presented a comprehensive framework enabUng assessment of the performance of joint 
spectrum sensing and MAC protocol operation for OSA networks. In the model we have proposed 
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we focused on the link layer throughput as the fundamental metric to assess performance. We have 
parameterized spectrum sensing architectures for energy detection based systems with collaborative 
measurements combining. We have proposed a novel spectrum sensing MAC denoted Truncated Time 
Division Multiple Access. We have also categorized multichannel MAC protocols for OSA networks 
based on their ability to buffer and switch existing SU connections on the arrival of a PU. Our analysis 
is supported by simulations which prove the accuracy of the obtained expressions. 

Some of the design guidelines that need to be noted are as follows. For spectrum sensing introducing 
TTDMA gives an improvement in obtained performance in compared to TDMA. Large networks, i.e. 
having many channels and users, benefit from clustering, while for small networks it is better to create 
small number of clusters such that sensing time is optimized. When considering MAC protocol design 
for OSA it is clear that more benefit comes from introducing SU connection buffering than channel 
switching, for those SU connections that have been preempted by PU. Interestingly, although intuition 
would suggest that MAC protocols that combine SU connection buffering and channel switching would 
outperform all other protocols, due to switching overhead this combination is usually inferior to protocols 
that involve only SU connection buffering. 

Our future task will be to investigate the delay experience by using any of OSA MAC protocols 
proposed. We plan to develop a comprehensive simulation software which will implement features not 
covered by our model, like queue per each SU. 
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